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Essentially all of our knowledge of the acceleration history of the Universe - including the ac¬ 
celeration itself - is predicated upon the validity of general relativity. Without recourse to this 
assumption, we use SNela to analyze the expansion history and find (i) very strong (5a) evidence 
for a period of acceleration, (ii) strong evidence that the acceleration has not been constant, (iii) 
evidence for an earlier period of deceleration and (iv) only weak evidence that the Universe has not 
been decelerating since z ~ 0.3. 


I. INTRODUCTION 


A still puzzling feature of our Universe is its accelerated expansion. Measurements of type la supernovae (SNela) 
provide direct evidence 0,110 and CMB anisotropy measurements provide important indirect evidence [3j. Within 
the context of general relativity, the acceleration can be explained if the present mix of matter and energy in the 
Universe includes a dark energy component that contributes about 2/3 of the critical density. Dark energy is defined 
by its large negative pressure, p ~ —p, and nearly spatially uniform distribution @j. 

While the energy of the quantum vacuum, mathematically equivalent to a cosmological constant, is the simplest 
explanation for the dark energy and cause of cosmic acceleration, all attempts to compute its numerical value have 
been unsuccessful, either leading to divergent results or numbers that are orders of magnitude too large; this is known 
as the cosmological constant problem 0 . Because it is possible that the energy of the quantum vacuum is zero or 
too small to explain cosmic acceleration, theoretical physicists have explored a plethora of possibilities for the dark 
energy, from a very light scalar field i i mum to the influence of extra dimensions E3 , and have even considered 
the possibility that there is no need for dark energy and that it arises due to heretofore neglected “ordinary” physics 
within general relativity (l3| . A fair summary of the present state of affairs is that there is little understanding of 
why the Universe is accelerating. 

Given that state of affairs, it is not unreasonable to think more broadly, even entertaining the possibility that 
the explanation may well have to do with gravity theory itself |14L 1151 llfil (l7l ll/ll . After all, we are confident that 
Einstein’s theory must be modified to make it compatible with quantum mechanics, and perhaps cosmic acceleration 
has something to do with that modification. However, relaxing the assumption of general relativity means that we 
must abandon the Friedmann equation, which relates the expansion rate to the matter and energy densities. We 
might then wonder what we really know about cosmic acceleration - or if the Universe is really accelerating since 
almost everything we have inferred about it depends upon the assumption of the existence of dark energy. 

To be specific, dark energy is often described as a smooth energy component with constant equation-of-state 
{w = p/p = const ~ —1), where quantum vacuum energy is the special case w = — 1 . The acceleration history of the 
Universe is then simply fixed by the matter content (UIm) and the value of w. It begins with a period of decelerated 
expansion and then crosses over to accelerated expansion at redshift 
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see Fig. ^ Both u>CDM with w = —1.05 ± 0.14 and Om = 0.29 ± 0.03 and ACDM with Om = 0.30 ± 0.05 are 
reasonable fits to the SNe data and a wealth of other cosmological data |j. (This in fact is the basic evidence for a 
period of accelerated expansion.) The existence of an early period of decelerated expansion is also supported by the 
success of the gravitational instability theory of structure formation and of big-bang nucleosynthesis. In particular, 
unless the equations governing the growth of small density perturbations are radically different, it is not possible to 
grow structure from small perturbations in an accelerating universe 0. In addition, it is difficult for a universe that 
accelerates during BBN to reproduce the observed light element abundances unless one artificially adjusts with the 
cosmic neutrino asymmetry Il9j . 

In this paper we relax the assumption of general relativity, but retain the weaker assumption of an isotropic and 
homogeneous space-time (which is supported by ample observational data) described by a metric theory of gravity 
(because there are few if any viable non-metric theory alternatives). In this case, SNela still provide a direct kinematic 
probe of the expansion, and we show that we can directly infer the following: (i) very strong (5cr) evidence for a period 





2 



FIG. 1: Different parametrizations of q(z). The solid line is the best-fit 2-epoch model (Eq. the long-dashed line is the 
best-fitting linear model; the short-dashed line is the best-fitting ACDM model (S2 m = 0.31). The kinematic models suggest a 
later transition to acceleration than ACDM does. 


of accelerated expansion; (ii) strong evidence that the acceleration has not been constant; (iii) evidence for an earlier 
period of deceleration; and (iv) only weak evidence that the Universe has not been decelerating since z ~ 0.3. 


II. FRIEDMANNLESS COSMOLOGY 
A. Kinematics 

With only the assumption that space-time can be described by a metric theory that is on the large isotropic and 
homogeneous, the Robertson-Walker metric still pertains, 

ds 2 = —dt 2 + R(t) 2 [dr 2 + r 2 dfi\ , (1) 

where R(t) is the usual cosmic scale factor and for simplicity spatial flatness has been assumed. Likewise, the cosmic 
scale factor at a given epoch is related in the usual way to the redshift that a free-streaming photon has suffered since 
that epoch: R(t) = 1/(1 + z). The epoch dependent expansion and deceleration rates are defined as usual and are 
related to one another by 


H(z) = 

q(z) = 


R 
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R =lff-i-l 


RH 2 dt 


H(z ) = Hq exp 


[1 + g(u)]dln(l + u) 


( 2 ) 

(3) 

(4) 


The relationship between detected energy flux T of a source at redshift 2 and its intrinsic luminosity C is also 
unaffected, 


T = 


C 


4t rd£ 


dL = (1 + z) f 
Jo 


du 


(5) 

( 6 ) 


where is the luminosity distance. Defining the distance modulus /i for a standard candle (or standardizeable candle, 
e.g., SNela) in the usual way, 


fi(z) = [ m B {z) - M b ] = 51og(di/Mpc) + 25 


( 7 ) 
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the observable, n{z), is related to the acceleration history, q(z), 

1 + z 


H(z) = 25 + 5 log 


H 0 • Mpc J o 


du exp I — / [1 + q(v)]dha.v 


( 8 ) 


Here Mg and ms{z) are respectively the absolute and apparent magnitudes of the source. Eq. [5] provides the 
fundamental kinematic relation that relates the deceleration history to SNela measurements. 


B. Simple Kinematic Models 


Even with the structure provided by the Friedmann equation, it has been a struggle to find a generic way to describe 
dark energy; most have settled upon using the equation-of-state w. Relaxing the assumption of the Friedmann equation 
exacerbates the problem. The kinematic parameters, H and/or q(z) have been used; and it has been suggested >2£j 
that jerk, j = — ( d/a)/H 3 , might be a good kinematic choice since the case j = 1 corresponds to a model that goes 
from matter-dominated behaviour to consmological constant dominated behaviour (of course, there is good evidence 
that our universe began from a radiation dominated phase which might have even been proceeded by a de Sitter 
phase). Operationally, H , q and j are all equivalent since each one can be expressed using integrals or derivatives of 
the others. 

Since q is related to the equation-of-state in general relativity, q = [fijvr + 3w(l — f2jvf)]/2, and there doesn’t seem to 
be strong motivation to go beyond q , we focus here on the acceleration history, q(z). We will consider three approaches 
to parameterizing q(z): (i) linear expansion, q(z) = qo+zdq/dz; (ii) piecewise constant; and (iii) a principal component 
analysis. For all these analyses we will use the gold set of SNela culled by Reiss et al. which contains 157 well studied 
SNela between z = 0.1 and z = 1.76 pTj . 

Reiss et al. analyzed their own gold set using the linear expansion and found that qo < 0 at the 99.2% level, that 
dq/dz > 0 at the 99.8% level and that the error ellipse for these two parameters is highly degenerate. Based upon this 
they concluded that a transition from recent acceleration to past deceleration occurred at a redshift of z t = 0.46 + 0.13. 
Here we add that the best-fit parameters give a 12% goodness-of-fit ( qo = —.74, dq/dz = 1.6, y 2 /DOF = 175/154), 
on par with the best ACDM model, which has an 11% goodness-of-fit (11 m = 0.31, y 2 /DOF = 177/155). The fact 
that this model finds g 0 < 0 gives us confidence that the Universe has accelerated recently. Furthermore, the positive 
slope is certainly an indication that q(z) was higher in the past. But can we trust that there was ever a transition? 
Unfortunately, this linear model always implies a transition redshift when its two parameters have opposite signs. 
Hence, it is unclear whether the transition is real or just an artifact. To illustrate, suppose that Riess et al. had found 
Zt = 3; since their farthest SNela has z = 1.76, it would have been clear that zt was implied not by the data but by 
the model itself. Since Zt actually does lie in the observed redshift range, evidence for a transition exists but remains 
unconvincing. 

An alternate parametrization is a piecewise constant acceleration with two distinct epochs )22H : 


q(z) 


q 0 for z < z t 
qi for z > z t , 


(9) 


where qo and qi correspond to the average values of q in their respective epochs. Fig. [5] shows our confidence contours 
in the qo~qi plane for various values of z t . The best fit model has z t = 0.08, qo = —2.0, q\ = 0.11, y 2 /DOF = 173/153 
and a 13% goodness-of-fit; however, the y 2 minimum in Zt is broad and the constraints on qo and q\ vary significantly 
with Zt- Moreover, in the limit of Zt <C 0.1 (z t 1), the model becomes essentially a one-zone model as there are no 
SNe to constrain qo ( q ±) and the error on qo (qi) blows up. Marginalizing over Zt with a prior that takes this into 
account, zt € [0.1,1], one obtains the final constraint contour in Fig. [3 Again there is strong evidence for recent 
acceleration (negative qo contains over 99.9% of the probability) and weaker evidence for past deceleration (positive 
qi contains 91% of the probability). The marginalized distribution for z t is very broad (see Fig. |3J) because qo and q\ 
are poorly constrained for low and high Zt, and meaningful constraints to Zt cannot be obtained. 

As shown in Fig. ^ the best-fitting kinematic models appear different from ACDM - in particular, they both find 
a lower redshift for the transition to acceleration. To address the significance of this, we generated 1000 mock SNela 
data sets, assuming ACDM with Dm = 0.31 and each with 157 magnitudes and errors taken from the gold set. We fit 
both kinematic models to the simulated data sets and tabulated the best-fit parameters. Fig. 0and Fig. Q]show that 
the kinematic parameters we expect to measure in a ACDM Universe are consistent with the constraints given by the 
gold set. If there is a difference between the actual acceleration history of the Universe and ACDM, these kinematic 
models and the present data do not have the statistical power to demonstrate it. 
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FIG. 2: Confidence contours in the qo — q\ plane for the two-epoch model. The 68.3%, 95.4% and 99.7% contours are shown 
for various values of zt■ The \ 2 function is broad in the zt direction and minimized for zt = 0.08. The final frame shows the 
contours after marginalization over zt £ [0.1,1]. 



FIG. 3: Constraints on the parameters of the 2-epoch model are compared to their expected measurements in a ACDM Universe. 
The solid lines are probability density functions obtained from the gold set and the histograms were generated from 1000 mock 
datasets as described. The qo (qi) function and histogram were marginalized over qi (qo) and Zt £ [0.1, 1[. The Zt function and 
histogram were summed over the full range of qo and q\. 
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FIG. 4: Constraints on the parameters of the linear model are compared to their expected measurements in a ACDM Universe. 
The solid lines are probability density functions obtained from the gold set and the histograms were generated from 1000 mock 
datasets as described. The transition redshift is not an independent parameter; it is given by Zt = —qo/(dq/dz). 
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FIG. 5: Maximum goodnesses-of-fit for three-epoch q(z) models that are not accelerating today. Setting q(z < zt) = 0 and 
q(z > z e ) = 0.5, the solid, long-dashed and short-dashed lines show 2 e = 0.5, 0.3 and 0.15 respectively. 

Both the 2-epoch and linear q(z) models provide direct evidence for a recent period of acceleration and weaker 
evidence for a past period of deceleration. However, as discussed, they each have their shortcomings. To further 
illustrate, we ask the provocative question, could the universe be decelerating today? First, as a practical matter, 
the Hubble flow only begins to dominate local peculiar motions for 0 > 0.01; this implies that any knowledge of 
the expansion dates back to at least O.OliZ,^ 1 ~ 0.1 Gyr - so as far as we know, the Universe could have stopped 
expanding 100 Myr ago! Next, consider the standard Taylor expansion of the luminosity distance, 

d L = Hq 1 z+ i(i - q 0 )H„ 1 z 2 + ... (10) 

Because the acceleration term is second order in z, one is further limited in detecting its effect by the accuracy of 
measuring dL- 10% uncertainty in Hq dL implies that objects at redshift greater than z ~ 0.1 are needed to measure 
qo with an error of order unity. Thus, direct evidence of acceleration/deceleration more recent than 1 Gyr ago will be 
difficult to obtain; this can be seen quantitatively in Fig. [3 where for z t <0.1 the error ellipse on qo blows up. 

The reason for asking whether or not the Universe is accelerating today is not just academic; there exist dark energy 
models that predict a history of alternating acceleration and deceleration [2J, |2J| . To address this issue we consider 
a three-epoch model, 


do 

for z < Zt 


dl 

for zt < z < z e 

(11) 

0.5 

for z > z e , 



motivated by the evidence for a recent epoch of acceleration and an earlier epoch of deceleration. The parameters 
qo and zt allow us to address the question of whether or not the Universe has been accelerating recently. To do so, 
we set qo = 0; Fig. [5] shows our goodness-of-fit as a function of z t after minimizing y 2 with respect to h and q-[. 
With an early transition at z e = 0.3, a long epoch of recent deceleration is consistent with the data at the 10% level. 
Larger values of qo will decrease the goodness-of-fit, but note that there may exist better-fitting parametrizations 
that raise it. In summary, the present SNela data cannot rule out the possibility that the Universe has actually been 
decelerating for the past 3 Gyr (i.e., since z = 0.3). 


C. Principal Component Analysis 

Here, we follow Huterer and Starkman |25l] and let the data itself guide our parametrization. By computing a Fisher 
matrix flci . Wf ], we can find accurately measured parameters and use them to form robust conclusions about cosmic 
acceleration. We start by dividing the redshift range of the SNela into N bins of width A z and write 

N 

q( z ) = 5>«(*) 

2=1 


( 12 ) 
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where the Ci(z) equal unity inside the *th bin and zero outside (NB: as N —> oo, arbitrary q(z) can be specified). 
From the gold set and Eq. 0 it is straightforward to numerically calculate the (N + 1 )x(JV + 1) Fisher matrix F for 
the /3i and the Hubble parameter h. For simplicity, we calculate F using f3i = 0 for all *; however, we have found 
that F is rather insensitive to the choice of fiducial parameters. After marginalizing F over h , we find a new basis of 
orthonormal functions et(z) which diagonalize F. The deceleration rate can be expanded in terms of these modes or 
“principal components”: 


N 

q(z) = ^2aiei(z). (13) 

i= 1 

The eigenvalues of E -1 estimate the uncertainties <T(a,;) in the coefficients aq. The advantage of the principle com¬ 
ponent basis is that the errors on the oq are uncorrelated. With A z = 0.05, the 6 modes of q{z) with the smallest 
errors are shown in Fig. 0 

We can use the principal components to reconstruct q(z), i.e. using Eq. 01] as our model, we can fit the aq to 
SNela data. Since F is insensitive to fiducial parameters, we may keep the er(aq) as error estimates. Note that if F 
were sensitive to fiducial parameters, we would have to ensure that the best-fit model was the same as the fiducial 
model before calculating the cr(aq); such a model can be found iteratively. Table [I] shows the results of parameter 
fits which use progressively more modes, starting with the most well-constained. The 1st and 2nd mode coefficients 
are respectively about 5-cr and 2-cr negative while higher modes are consistent with zero due to increasingly large 
variances. Increasing the number of modes slightly alters the fit; however, since the data is insensitive to higher modes 
and since the coefficients are uncorrelated, we find it unlikely that keeping more modes will change the coefficients in 
Table □ by 1 -a. 


# modes 

x 2 

h 

Qi 

«2 

a 3 

Q!4 

1 

177.7 

.644 

-1.10 

- 

- 

- 

2 

173.2 

.653 

-1.08 

-1.61 

- 

- 

3 

172.6 

.664 

-1.08 

-1.71 

-1.10 

- 

4 

172.6 

.662 

-1.08 

-1.56 

-1.30 

.61 

cr{a n ) 



0.22 

0.72 

1.52 

2.54 


TABLE I: Reconstruction of q(z) by Eq. II .'fl nsing well-constrained principal components. The top part of the table shows the 
results of parameter fits as progressively more modes are added. The estimated errors on the parameters follow directly from 
the Fisher matrix used to derive the modes. We have marginalized over h. 


It is sensible to reconstruct q(z) using only the well-measured modes J2f|. Doing so gives us a coarse yet accurate 
picture of the deceleration parameter in its “mode space,” but it also introduces a bias: the reconstruction will lack 
features of the true q(z) which are poorly probed by SNela. Consider for instance the two-mode reconstruction 
illustrated in Fig. 0 it approaches zero at high redshift, but it is clear that we cannot trust this feature. It is also 
clear that the reconstruction will look different if we include more modes. So which features of the reconstruction can 
we trust? 

First, we can trust that the Universe has accelerated. To see this, we note that the 1st mode is completely positive. 
Since the modes are orthonormal, the coefficients oq are given by 

a» = J dz q(z)e.i(z). (14) 

It follows that ai is negative only if q(z) < 0 for some z. Our estimate that oq is 5-cr negative therefore implies a 
5-<j detection of cosmic acceleration at some redshift. The shape of ei(z) implies that this redshift is likely to be low. 
This conclusion is not model-dependent since we are free to decompose q(z) using any complete basis we like. One 
could make a similar argument using e.g. a polynomial or Fourier expansion, but the principle component expansion 
is advantaged since the uncertainty in aq is small and non-degenerate with other parameters. 

Next, we can also trust that q(z) is not constant. Fitting a constant q model to the gold set yields q = —0.29 with 
a 6% goodness-of-fit, but fitting only to modes that have been convincingly measured shows that constant q really 
fails. Using the reconstructions in Table|l| we fit a constant q model directly to the mode coefficients. The best fit to 
these reconstructions is also q = —0.29; it is essentially determined by the first two modes, which have the smallest 
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FIG. 6: The six most well-constrained principal components of q{z). They are normalized so that ei(0) > 0 and J ei(z) 2 dz= 1. 
The errors on the mode amplitudes, o-(cti), are estimated from the Fisher matrix. 



FIG. 7: Reconstruction of q(z) using its 2 most well-constrained principle components (Eq. 11 31 . The short-dashed line is the 
best 2 mode fit to SNela data; it has a 14% goodness-of-fit. The gray band shows the error corresponding to 1 -a uncertainties 
in both ai and « 2 . The solid line is q{z) for ACDM with £Im = 0.3. As in Fig. Q the kinematic model suggests a later 
transition to acceleration than ACDM does. 


























Mode number n 


FIG. 8: Mode space representation of q{z). The solid line indicates the a n corresponding to q(z) = —0.29. The points show 
the values of a n obtained from a 4-mode fit, and the error bars are taken from the Fisher matrix (see TableCJl. 


errors. We obtain the following goodnesses-of-fit: 


P(X 2 > Xmin) 


0.1% first 2 modes 
0.4% first 3 modes 
1% first 4 modes 


(15) 


The main reason for the poor fit is easily seen in Fig. 0 if q were constant, the first two coefficients should have 
roughly opposite values, but they don’t. Higher modes will tend to improve the fit by adding degrees of freedom 
without significantly increasing y. 


III. CONCLUSIONS 

Understanding dark energy may ultimately require us to abandon or change the familiar equations of general 
relativity. Hence it is crucial that we keep track of which observations and conclusions depend on the assumption 
of general relativity and which are more robust. Without the Friedmann equation, it is still possible to describe 
the Universe with a Robertson-Walker metric and infer the expansion history via the dimming and cosmological 
redshifting of SNela. We used this fact to learn about cosmic acceleration kinematically. 

In particular, we have shown that the gold set of SNela from Riess et al. provides the following information: 

• Very strong evidence that the Universe once accelerated — Though various parametrizations can imply 
this, the strongest evidence comes from measuring the best-constrained principal component of q(z). The 
acceleration is likely to have been relatively recent in cosmic history. 

• Strong evidence that q was higher in the past — The constant q model has a poor goodness-of-fit, especially 
when fit in mode space using only well-measured modes. In addition, the linear q{z) model indicates that the 
average dq/dz is positive. 

• Weak evidence that the Universe once decelerated — Both the two-epoch and two-mode models of q{z) 
slightly favor past deceleration. The linear model also predicts a transition to positive q. Alas, it is difficult to 
show that this is not a model-dependent feature. 

• Little or no evidence that the Universe is presently accelerating — It is manifestly difficult to constrain 
q(z < 0.1) with SNela. Using a three-epoch model, we further showed that a long period (e.g., z = 0 ^ 0.3) of 
recent deceleration is consistent with data. 

• ACDM and wCDM are acceptable fits to SNela data. While there are tantilizing indications that the 
actual acceleration history may be different, none of the kinematic models studied here have revealed robust 
features about cosmic acceleration that differ from ACDM or wCDM. Moreover, none of the kinematic models 
fit the data significantly better. 
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